Out-of-equilibrium critical dynamics at surfaces: Cluster dissolution and 

non-algebraic correlations 
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We study nonequilibrium dynamical properties at a free surface after the system is quenched from 
the high-temperature phase into the critical point. We show that if the spatial surface correlations 
decay sufficiently rapidly the surface magnetization and/or the surface manifold autocorrelations 
has a qualitatively different universal short time behavior than the same quantities in the bulk. 
At a free surface cluster dissolution may take place instead of domain growth yielding stationary 
dynamical correlations that decay in a stretched exponential form. This phenomenon takes place in 
the three-dimensional Ising model and should be observable in real ferromagnets. 



Dynamical properties of systems which are quenched 
from the high-temperature phase to or below the critical 
point are of recent interest, both experimentally and the- 
oretically (see RefiSi* for reviews). In the bulk of the 
system, which is generally studied in this context, an or- 
dering phenomenon takes place, which is manifested by 
the existence of a characteristic length-scale, £(i), grow- 
ing in time as £(t) ~ t 1 ^, where z is the well known 
dynamical exponent. If the quench is performed to the 
critical temperature new (nonequilibrium) dynamical ex- 
ponents have to be defined. The short time dependence 
of the magnetization starting from a small finite initial 
value, nii > is given by m(t) ~ rriit 0b , where the initial slip 
exponent is expressed as* 9b — (xi — Xb)/z . Here Xb is 
the static scaling dimension of the bulk magnetization, 
whereas Xi is a new nonequilibrium exponent, the scaling 
dimension of the initial magnetization. Since generally 
due to fluctuations Xi > Xb there is a universal small- 
time increase of the order. The initial disordered state 
has a long-lasting effect on the autocorrelation function, 
too, measured at time t, after a waiting time s < t. In 
the limit t>s the autocorrelation function behaves as* 
C(t,s) ~ t~ Xb / z , where the autocorrelation exponent in 
a d-dimensional system is given by A& = d — Xi + Xb- 

Real materials are bounded by surfaces and many 
nonequilibrium processes (thermalization, transfer of 
heat, etc.) are very intensive at the surface. At a free 
boundary layer the order is weaker than in the bulk due 
to missing bonds and surface correlations involve a new 
scaling dimension, x\ > Xb- Due to the inhomogeneity 
at the surface the spatial and dynamical correlations are 
anisotropic: critical correlations between a surface point 
and one in the bulk has an exponent x± + Xb and only 
correlations between two surface points involve the expo- 
nent 2xy&l. 

The phenomena of nonequilibrium dynamics have also 
been studied at surfaces by numerical and by scaling and 
field-theoretical (FT) methods*!*. It is found that the 
scaling dimension of the initial magnetization, Xi , is the 
same as in the bulk, but the local magnetization operator 
has its surface dimension, x±, instead of Xb- Therefore 



the exponent of the relaxation process, which involves 
one magnetization operator, is expected to change as 
9i = 8b — (x\ —Xb)/z, whereas the autocorrelation of two 
operators will have an exponent** of Ai = A/, + 2(xi — Xb). 
The FT results have been checked numerically, but the 
simulations were restricted to the two-dimensional (2d) 
and 3d Ising models. 

In this Letter we are going to revisit the problem of 
nonequilibrium critical dynamics at surfaces. We show 
that during domain growth (DG) the correlated sites in 
the domain are either percolating (to which we refer as 
conventional DG dynamics) or they form isolated finite 
clusters which gradually shrink as time goes on (in which 
case we use the term cluster dissolution (CD)), depending 
on the critical exponents of the given system. In the CD 
regime relaxation and short time autocorrelations have 
unconventional, stretched exponential time dependence, 
which are shown by scaling theory and checked by nu- 
merical calculations. This process takes place in the 3d 
Ising model and should be experimentally observable in 
real magnets. Surface manifold autocorrelations and per- 
sistence have also been studied. 

In a phenomenological theory the number of correlated 
points within a domain of volume V ~ £ d can be es- 
timated as the result of two competitive processes, i) 
Initial nonequilibrium growth creates correlations among 
Vi ~ £, Xi sites, ii) which are reduced by equilibrium crit- 
ical fluctuations by a factor of £~ x > . (Here we use the 
notation xi for the local exponent, which is Xb in the 
bulk and X\ at the surface.) If X{ > xi, the correlated 
points in the domain form a percolating cluster, the mass 
of which scales as: M ~ \ n this DG regime the 

magnetization at time t is given by m(i) ~ rriiM with 
the initial magnetization m t . This formula recovers the 
known FT results, both in the bulk and at the surface. 
For Xi < xi, however, the correlated points in the do- 
main form only finite clusters and in this regime the CD 
dynamics applies. 

The autocorrelation function in the long time limit, 
s > 1, i - s » s, is invariant under the scale trans- 
formation, C(t,s) = b~ n C(tb~ z , sb~ z ), when lengths are 
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rescaled as L — » L/b, b > 1. Here 77 = 2a;; is the decay 
exponent of the local magnetization. With the choice 
b = £ s = s x l z we arrive at the result: 



C(t,s) = c(s)f(t/s), c(*)~£ 



t - s > s 



(1) 



where the pre-factor, c(s), can be interpreted as the prod- 
uct of the critical magnetizations in two cells of size £ s . 
In this interpretation the domain at time t 3> s, with typ- 
ical size £, is divided into a large number of cells. In the 
short time limit, t — s < s, however, when the domain 
contains only one cell, the spatial correlations between 
two points of reference are in the DG regime of local-bulk 
type, leading to the scaling form of the autocorrelation: 

C(t,s)=c 1 (s)f 1 [(t-s)/s}, Cl (s)~C m , t-a<s, 

. (2) 

with rji = Xb + x\ . Thus in an inhomogeneous environ- 
ment the autocorrelation function can be formally writ- 
ten as: 



C(t,s)=c(s)f(t/s) 



c(t) 
ci(t) 



Ci(M), t-s»s, (3) 



where C\{t 1 s) is the autocorrelation in a homogeneous 
environment with local exponent xi. In the DG regime 
C\(t, s) follows from the observation that out of the vol- 
ume of the domain there are M sites which are correlated 
with the initial one, leading to 



d(t,s) ~ M/V ~ C (d ~ Xl+Xl) ■ 



(4) 



For a bulk site, with x\ = Xf, and c{t) — ci(t) we im- 
mediately obtain C(t,s) ~ f-(d-xi+x b )/ z m acc0 rdance 
with FT calculations. Similarly, for a surface point with 
xi = x\ and c(t)/c\{t) ~ t^ Xb ~ Xl ^ z we recover the FT 
result. 

Let us now turn to the CD regime, when i) the cells 
contain only finite clusters and ii) in the short time limit 
the two points of reference are within the same cell. In 
a given cell the probability to find a large cluster of lin- 
ear size, L, is exponentially small: P(L) ~ exp(— aL d ), 
but the relaxation time, t r , related to this cluster is 
large and can be estimated as follows. Correlations 
within the cluster are destroyed if a domain wall of 
W ~ L d ~ x sites are created within it. The number of 
still correlated sites in the wall at time t is estimated as 
~ Wt^ Xi ~ x '^ z , from which the value of the relaxation 
time, t r ~ L z i d - 1 )/( x i- x i) ^ follows. Now the dominant 
behavior of the autocorrelation function is obtained as 
an average over the clusters and is given by: 



(5) 



C(t) ~ / P(t r )exp(-T/t r )dt r ~ exp(-CV K ), 

_ (xj - Xj)d 
K ~ z(d-l) ' 



in terms of r = t — s < s. If it happens that k > 1, then 
the decay is pure exponential. Note that the autocor- 
relation function is stationary as long as the two points 



of reference are within the same cell. In the long time 
limit, t>s, however, when the two points are in different 
cells, one can use a coarse-grained picture in which the 
DG theory, with the results given below Eq.(0}, becomes 
applicable. 

We can also estimate the short time relaxational be- 
havior of the initial magnetization, rrii, by noting that the 
dominant process in this case is the same as for the auto- 
correlation function. From this we obtain a fast stretched 
exponential short time decay of: 



m{t) ~ rrii exp(—C m t K ), t < t t 



(6) 



with the universal exponent, k, as given in Eq.lQ and the 
decay is pure exponential if K > 1. At ti ~ |lnmi| _1 / K 
the relaxation crosses over to a power-law behavior. 

The above scaling results are checked by numerical 
simulations of a 2d Ising model in which the strength of 
the coupling at a distance y from the free surface deviates 
asymptotically from its bulk value by an amount of A/y. 
In this Hilhorst-van Leeuwen (HvL) modelifi, the sur- 
face magnetization scaling dimension is known exactljsii: 
x\ — (1 — A/A c )/2, for A < A c and A c is a constant. The 
presence of the surface defect does not modify the value 
of the bulk exponents, which are known to be: Xb = 1/8, 
Xi = .53 and z = 2.17. 

In the actual simulation we considered L x L finite 
systems up to L = 300 with periodic boundary condi- 
tion (b.c.) in the direction of translation invariance and 
free b.c. in the perpendicular, non-homogeneous direc- 
tion and made measurements up to several hundred time 
steps. The calculated short time behavior of the surface 
magnetization is shown in Fig. la. It has indeed different 
characteristics for xi > x\ and for Xi < x\, respectively. 
In the former case, the initial slip behavior sets in after a 
short time period, whereas in the latter (CD) region there 
is a fast non-algebraic short time decrease of the magne- 
tization (grey lines) before the much slower power-law 
decay starts. In the fast short time region the decay cor- 
rectly fits the predicted universal stretched exponential 
behavior of Eq. ©, as shown in the inset of Fig. la. 

In our next calculation we compared the short time 
behavior of the surface autocorrelation function in the 
3d and in the 2d Ising models. For the 3d Ising model 
with Xb = .516, x\ = 1.26 (ordinary transition), Xi — 
.74 and z — 2.04 we are in the CD regime, since Xi < 
x\, whereas the (pure) semi-infinite 2d Ising model is in 
the DG regime. For the 3d model the simulations were 
performed on a cube of linear size up to L = 60 with free 
b.c. in one direction and periodic b.c.-s in the other two 
directions. The numerical results for 3d (Fig. 2a) show a 
stationary regime for t— s < s in which the decay fits very 
well with the theoretical prediction in Eq.© (see inset). 
In contrary for 2d (Fig. 2a) there is no stationary regime, 
even not for very short time differences (see inset). 

We have also considered nonequilibrium dynamics re- 
lated to a d! dimensional manifold, in particular in the 
case where the manifold is the surface layer of a semi- 
infinite system, thus d' = d — 1. The manifold autocor- 
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FIG. 1: Short time behavior of (a) the surface magnetization 
and (b) the surface manifold autocorrelation function G(t, s) 
at the waiting time t = s in the HvL model for various val- 
ues of the surface magnetization scaling dimension x\. The 
broken straight lines indicate the predicted long time asymp- 
totics. Here and in the following, the data, which have been 
obtained by averaging over typically 30000 different runs, are 
free of finite size effects. The inset in Fig. la shows the ex- 
pected stretched exponential behavior for x\ = 1 
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FIG. 2: Surface autocorrelation function of (a) the 3d and (b) 
the 2d semi-infinite Ising models as function of t—s. As shown 
in the insets, a stationary stretched exponential behavior is 
observed in 3d with the exponent k in Eq.(0, whereas in 2d 
no stationary short time regime is found. 



relation function, G(t,s), at the waiting time t = s has 
a non-trivial behavior: G(t,t) ~ ^ d _x i _a! i>) ) provided 
d'f = d! — xi — x\j > 0. Thus in this case the initially cor- 
related points in the manifold form a percolating cluster 
inside the cell, with a fractal dimension of d'j . This sit- 
uation corresponds to the DG process in this problem. 
If, however d! — xi — x\> < 0, the manifold autocorrela- 



tions G(t, t) approach a finite, time independent value, 
thus the initially correlated sites within the growing cell 
are finite, isolated clusters, which are then the subject of 
dissolution during the CD process. These predictions are 
successfully verified by numerical calculations on the 2d 
HvL model for different values of x\ as shown in Fig. lb. 

For t > s the considerations about G(t, s) are simple 
in the DG region, i.e. when d! — xi — Xb > and the 
result can be obtained by integrating the related formu- 
lae for a single spin. For example we get from Eq. for 
the asymptotic behavior of the manifold autocorrelation 
function in the long time limit: G(t,s) — g(s)f'(t/s), 
where g(s) ~ s( d ~ 2x i)/ z . Thus the waiting time depen- 
dence is a subject of change as time goes on, and the 
scaling combination: s( x '~ Xb '' z G(t, s)/G(s, s) is just the 
function of t/s, as illustrated in Fig. 3 for the 2d HvL 
model. The time dependence of the manifold autocor- 
relation function in the long time limit is obtained from 
Eq.© as G(i,s) ~ t- x '/ z . Here A' = d-d' -x l -x b + 2x ll 
from which both the bulk (x; = Xb) and the surface 
(xi = x\) results follow. 
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FIG. 3: Rescaled surface manifold autocorrelation as function 
of t/s in the DG regime of the HvL model. For illustrative 
reasons the data points for x\ = 1/2 resp. x\ = 1/4 have been 
shifted horizontally by 2 resp. 4 units. The expected stretched 
(pure) exponential behavior in the CD regime is shown in the 
inset. 

In the CD region, when d! — xi — Xb < 0, G(s, s) has a 
finite value, thus in the cell of size £ s there are only finite 
clusters. In the short time limit, t—s < s, as for the single 
spin case the two points of reference are within the same 
cell and the autocorrelations are governed by the CD 
dynamics. In this case the driving force of the dissolution 
is the inhomogeneity effect as local-bulk correlations go 
in time over to local-local correlations, as seen in Eqs.QJ 
and (J2J) for a single spin. In a time period, r = t — s < s, a 



cluster shrinks by a factor of ~ r 



(x b ~xi)/z 



Now using the 
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same reasoning as for the single spin problem we arrive 
to the same conclusion, just Xi in the single site result 
has to be replaced by x b for a manifold. Thus the surface 
manifold autocorrelation function in the short time limit, 
t — s < s, is stationary and described by a stretched 
exponential form: G(t, s) ~ exp(— C'(t — s) K ), with k' = 
(xi — Xb)d/(z(d — 1)), unless k! > 1, when the decay is a 
pure exponential. Again, in the long time limit, t — s 3> s, 
when the two points of reference are in different cells in 
a coarse-grained picture we recover the conventional DG 
results. Numerical calculations for the 2d HvL model are 
in complete agreement with the theoretical predictions. 
As seen in the inset of Fig. 3 the decay is stationary 
and for x\ — 1 it is stretched exponential, whereas for 
x\ = 5/4 it is pure exponential. 




FIG. 4: Time dependence of the surface manifold persistence. 
The main figure shows the stretched (pure) exponential be- 
havior in the regime d! — xi — x b < 0, whereas the power-law 
behavior for d! — xi — xt > is displayed in the inset. 

We close our Letter by presenting results concerning 
the surface manifold persistence, which completes recent 
studies by Majumdar and Brayi^ about persistence of 
manifolds which are embedded into the bulk. Since per- 



sistence and (the short time behavior of) manifold auto- 
correlations are intimately related, the persistence also 
shows qualitatively different behavior for dl — xi — Xb > 
and d' —xi~x b < 0, respectively. In the first (DG) regime 
the persistence is in a power-law form: P pr {t) ~ i~ epr , 
where Q pr is some non-trivial function of d'^ = d! — xi— x b . 
This result follows from a general theorem by Newell and 
Rosenblatt. On the other hand in the CD regime with 
d! — xi — Xb < the persistence has a faster, stretched ex- 
ponential or pure exponential behavior, as the autocorre- 
lation function: thus P pr (t) ~ exp(— ait K ), if < k' < 1, 
otherwise P pr (t) ~ exp(— a^b), if k' > 1. We note that 
for a bulk manifold the persistence shows also these three 
different types of functional forms, however the border 
between the different regions as well as the actual value 
of the exponent in the stretched exponential are different 
in the two problems^. 

Numerical results in Fig. 4 for the 2d HvL model are in 
complete agreement with the predicted stretched (pure) 
exponential decay in the regime d! — Xi — Xb < 0. Results 
in the d! — xi — Xb > regime are presented in the inset to 
Fig. 4. The long-time power-law dependence is clearly 
seen for the different models, which include: i) the surface 
manifold of the 3d Ising model, ii) the surface manifold of 
the 2d Ising model {x\ = 1/2) and iii) the case x\ — Xb — 
1/8 where we expect the persistence of a bulk manifold 
in the 2d Ising model. 

In summary, we have presented a scaling theory de- 
scribing nonequilibrium critical dynamics at free surfaces. 
Depending on the values of the involved scaling dimen- 
sions, domain growth or cluster dissolution is observed. 
The latter phenomenon, which has not been observed up 
to now, leads to unconventional behaviors of dynamical 
correlation functions. The theoretical predictions have 
been verified in extensive Monte Carlo simulations. 
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